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m ABSTRACT 

o 

According to a standard ohmic perspective, the injection of accelerated elec- 
trons into the flaring region violates local charge equilibrium and therefore, in 
response, return currents are driven by an electric field to equilibrate such charge 
violation. In this framework, the energy loss rate associated to these local cur- 
rents has an ohmic nature and significantly shortens the acceleration electron 
path. In the present paper we adopt a different viewpoint and, specifically, we 
study the impact of the background drift velocity on the energy loss rate of ac- 
celerated electrons in solar flares. We first utilize the Rutherford cross-section 
to derive the formula of the energy loss rate when the collisional target has a 
finite temperature and the background instantaneously and coherently moves up 



Oh 



O 

to equilibrate the electron injection. We then use the continuity equation for 
i_i electrons and imaging spectroscopy data provided by RHESSI to validate this 

,__i model. Specifically, we show that this new formula for the energy loss rate pro- 

vides a better fit of the experimental data with respect to the model based on 
^ the effects of standard ohmic return currents. 
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Introduction 



In arc-shaped solar flares, huge amounts of electrons are accelerated due to magnetic 
reconnection. In the standard picture (Fletcher et al.|2011 ) accelerated electrons are injected 



at the top of the loop-shaped flare, move down along magnetic field lines and, during their 
motion, loose energy because of various mechanisms. Coulomb collisions with ambient parti- 
cles represent the most relevant mechanism explaining energy losses for accelerated electrons 
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in loop-shaped flares. In this framework, the predominant energy loss process involves the 
interaction between the injected electrons and the electrons in the flaring target. Specifi- 
cally, under the assumption of a cold target, the energy loss rate formula depends as 1/E 



on the electron energy E (Emslie 1978) while in a warm target a more complicated formula 



accounts for the temperature T of the ambient medium ( Spitzer|1962 Longmire|1963 ). How- 
ever, other important effects may impact on this process, such as wave-particles interactions 



(Hoyng & Melrose 1977, Hannah & Kontar 2011) and return currents. This paper focuses 



on the role of return currents in the energy loss mechanism during flares and specifically 
utilizes measurements observed by the Reuven Ramaty High Energy Solar Spectroscopic Im- 
ager (RHESSI) to empirically study the mechanism of energy loss in the presence of return 
currents. 



In their pioneering paper, Knight & Sturrock (1977) analyzed the effect of return cur 



rents on accelerated electron motion by means of electrostatic fields and using Vlasov equa- 
tion. The combination of reverse current electrostatic field with collisions has been discussed 



in 



Emslie (1980), where the impact of return currents on hard X-ray emission has been 
studied as well. Bounds for the size of the electrostatic term due to unstable electron-ion 



drifts are found in Emslie (1981). In Brown & Bingham (1984), another approach using 



electrostatic modeling is presented. The influence of inductive fields and their interplay with 
the electrostatic one have been thoroughly discussed and compared by |Larosa fc Emslie 



(1989) and by van den Oord (1990). Starting form observations, McClymont & Canfield 



(1986) argued that the areas of injection are considerably small and thus very large return 



currents are expected. More recently Alexander & Daou (2007) used RHESSI observations 



to estimate the injection area, confirming the effects on hard X-ray emission predicted in 



Emslie (1980). A complete treatment of the many possible kinetic energy loss mechanisms, 



including ohmic losses, has been more recently presented in Zharkova & Gordovskyy (2005). 



All previous papers are based on an ohmic viewpoint according to which the injection 
of a large number of accelerated electrons into the loop violates the local charge equilibrium 
and, in response to this violation, local currents are established by flare background electrons. 
In this picture, the generated return currents are driven by an electric field and energy losses 
of ohmic nature shorten the accelerated electron path before thermalization. However a 
different model describing the restoration of charge equilibrium is possible, in which the 
production of return currents occurs instantaneously, the background particles moves up 
coherently along the field line, while the energy loss rate due to Coulomb collisions is modified 
in a very peculiar manner: since the single accelerated electron appears more energetic in the 
rest frame of the background motion and since the energy loss rate decays with energy, we 
expect the resulting Coulomb collision energy loss rate to be lower than the one computed 
for background particles with vanishing drift velocity. 



The aim of the present paper is to investigate the impact of this coherent background 
motion on the energy loss rate for accelerated electrons and to compare it with the compet- 
ing influence of ohmic forces induced by return currents. This comparison will be performed 
by following the same empirical approach adopted by Torre et al. (2012): we will (1) 



use 



hard X-ray imaging spectroscopy data observed by RHESSI to reconstruct mean electron 



flux images of extended sources (Piana et al. 2007); and (2) apply the electron continuity 
equation to select the model for energy loss rate that best fits the empirical electron maps. 



The plan of the paper is as follows. In the next section we shall derive the energy loss 
rate for Coulomb collisions when a coherent motion is added to the background Maxwellian 
distribution. In the third section we shall describe the procedure we have adopted to com- 
pare the various energy loss rate models. The comparison of such models with RHESSI 
observations is presented in the fourth section and finally some conclusions are drawn in the 
last section. 



2. A general model for the Coulomb collisional energy loss rate 



The main ingredients at the basis of our model are two formulas obtained by Rutherford 



(1911) and Butler & Buckingham (1962), respectively. In the first one the differential cross 



section, which describes the probability of a single collision of a particle of charge ze- and 
mass m travelling with a speed v against a target particle of mass M and charge Ze_ moving 
at velocity w, is 
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This expression is non zero only for scattering angles 9 bigger than a fixed small angle 9q 
which is related to the Debye screening length. 

The second formula describes the collisional loss rate for a particle traveling in a plasma, 
whose particle velocities are distributed as /(w), and is given by 
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In order to discuss the role of a coherent motion in the energy loss rate, we shall follow 
the derivation of pi) given by Butler & Buckingham ( 1962 ) and modify the form of the target 



distribution. We shall here assume that the background particles are in a stationary solution 
of the Boltzmann equation, i.e. that they are distributed according to the Maxwell law but 
with a coherent motion induced by a drift velocity wo, and a thermal velocity wt such that 
Mw T /2 = ksT ', where T is the temperature of the distribution and ks is the Boltzmann 
constant. Under these hypotheses we have 
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where n is the target electron density. Integrating over the target velocities and expressing 
the result in terms of v' = v — w we get 
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where, K = 2ir(Zz) 2 e 4 _A and erf is the standard error function. It is possible to point out 
two terms in equation (J5J) characterized by two different physical meanings. The first term 
corresponds to the pure thermal energy loss rate computed in a frame where the coherent 
velocity of the background vanishes (w = 0). The latter one is proportional to v' ■ w 
and cannot be directly inferred from the pure thermal energy loss rate. Equation (J5J) can 
be simplified by assuming that wo and v are anti-parallel and that the electron-electron 
collisions give the most important contribution to the energy loss. Therefore, if 
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When wo vanishes the previous expression reduces to the standard energy loss rate formula 



that can be found in Spitzer (1962). As previously noticed in the general case, also with the 



further assumption stated above, the energy loss rate is formed by two contributions. The 
one which cannot be directly inferred from the thermal energy loss rate formula is relevant 
also in the zero temperature limit (% — > 0) where we obtain 
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Fig. 1. — The energy loss rate for the electron-electron cross section with thermal energy 
fixed at 1.6 keV for various background velocities Wq. In the plot we have used as labels the 
kinetic energy E associated to the corresponding drift velocity Wq. The red line indicates 
the case when E = 0. 



Figure [I] provides a qualitative analysis of the energy loss rate ^ for different values of the 
kinetic energy E associated to the drift velocity Wq. From this analysis it clearly follows that, 
when the intensity of the background motion increases, the energy loss rate of the injected 
electrons decreases. Further, all different forms of the energy loss rate for all different values 
of Eq assume an asymptotic behavior 1/E at high electron energies. 



3. Selection of the energy loss model 

We now validate the model in equation ([7]) with electron flux maps reconstructed from 



RHESSI data, following the approach introduced by Torre et al. (2012) . 



In the standard flare picture, electrons move from the injection region towards foot- 
points along magnetic field lines, loosing energy along the paths. If F(E, s) (electrons cm -2 
s _1 keV -1 ) is the electron flux, differential in energy along the direction s and N(s) is 



the column depth, Emslie, Barrett, & Brown (2001) proved that g(s; E) = N(s)F(E,s) 
(electrons cm" 4 s _1 keV -1 ) satisfies the continuity equation 
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where the first sign is positive for electrons moving towards larger s and negative otherwise. 
In this equation, g(s; E) corresponds to the mean electron flux maps that can be recon- 
structed from RHESSI visibilities ( Piana et al.|2007 ) while S(s, E) is the source term coding 
the information on the injection region. The continuity equation can be interpreted to obtain 
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is an empirical quantity determined from the electron maps. In (pi) and (J9J)-(10) the source 
term S(s, E) describes the injection of electrons in the flare region and thus it is an energy 



gain term. Following Guo et al. (2012) here we assume 
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where h s (electrons cm -5 keV -1 s _1 ) is the source amplitude averaged along the line of sight, 
5 is the spectral index of the injected electrons, L represents the length of the injected region 
and E s is fixed equal to 10 keV. Since here we are interested in the energy domain, and since 



we shall average over different s as in Torre et al. (2012), the choice of a box shaped injection 



as in (11) is not really restrictive. 



The total energy loss rate 



\ ds J tot 
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in (J9J) represents the core of the present analysis. In the following we will consider four 
possible situations: 

• Model 1: hot target with charge equilibrium. This is the case pictured by (J7l), in which 
the collisional target has a finite temperature and the background instantaneously 
reacts to the electron injection by means of a coherent motion. 



Model 2: cold target with charge equilibrium. This is what happens when in ([7]) one 
fixes E T = \mw% = 0. 



Model 3: ohmic losses. In this case there is no drift velocity in the collisional term 
([7]) and the energy loss rate has the form 
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where the second term on the r.h.s. describing the ohmic losses is given by 
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Here £ is the electric field which drives the return currents, rj is the resistivity, n is the 
background density and wq is again the background velocity. The first term at the r.h.s. 



is the standard energy loss formula that can be found in Spitzer (1962), corresponding 



to the case when no background velocity is present and the target temperature is finite. 



Model 4: hot target without return currents. This is the model described by Spitzer 



(1962) and already discussed by Torre et al. (2012). 



For all four models the number of free parameters to fit against the empirical R(E, s) is 
always two, and in all cases these parameters are the target density n and the averaged 



source amplitude h s in (11). In fact, the target temperature can be fitted by using spatially 



integrated spectroscopy. This same spectroscopy and the charge equilibrium constraint allow 
fixing the drift velocity Wq. More specifically, assuming that the injected electron flux 
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is equal to the flux nN J ^E associated to the return currents leads to 
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Since the emission measure EM can be inferred from spectroscopy and the flare area A from 



the reconstructed electron maps, equation (15) represents a constraint for Eq and therefore 
for Wq- 



4. Models versus observations 

In this section, we present the results of the analysis stated above for various time 
bins of three different flares, namely those occurred on April 15 2002, April 16 2002 and 
May 21 2004. The time intervals, and the relevant physical parameters are described in 
the following Table 1. The temperature T, the spectral index 7 = 5 — ijj and EM are 
obtained from photon spectral fitting. The flare area A is computed from the reconstructed 
electron maps. Specifically, we used the map at 14 keV as reference image and summed up 
the areas of pixels with intensity higher than 10% of the maximum intensity. The spatial 
region As = [s min , s max ] over which the average of R(E,s) is computed is determined as in 



Torre et al. (2012) (see Figure 4 of that paper). Specifically, the electron flux images are 
considered for all energies and for each row of each image we took the pixel with maximum 
intensity. We thus drew a path in each image approximating a field line. In this path we 
fixed at s = the pixel with maximum intensity and assumed that moving toward the right 
foot-point increases s while moving toward the left foot-point decreases s. 

The values of h s and n are determined as described in the previous section, i.e. by 
fitting the empirical R(E, s) values deduced from the mean electron flux spectral images and 
averaged along As, against the four models of the energy loss rate. These best-fitted values 
for h s and n in the three events under analysis are given in Table 2 where we also provide 
the corresponding values of the kinetic energy associated to the return currents' velocity. 
Furthermore, Figure |2J contains the empirical values of R(E) used for the fitting and shows 
the best-fit curves corresponding to the four theoretical models. 

From these results we notice first that for the April 16 2002 event (in both considered 
time intervals), the x 2 values corresponding to the four models are very similar. Correspond- 
ingly, the values of the kinetic energy Eq associated to the drift velocity w$ obtained for this 
event are considerably smaller than for the other two events. This is probably a consequence 
of the fact that this event in these time intervals is in a late phase (coherently, the panels 
in Figure ^ corresponding to this event show decreasing values for R(E,s)). For the other 
data sets, Model 1 systematically provides smaller \ 2 values. Knowing the emission measure 
obtained from photon spectral fitting, and that EM = Vn 2 , where V is the flare volume, 
we can validate the obtained density values. Assuming that the flare has the shape of a 
tube, we can estimate V for every time intervals starting from the electron flux maps as 
previously done for the area A. We obtain that the emission measure obtained for Model 1 



1 This relation holds exactly in the case of the Kramers bremsstrahlung cross-section, while for more 
general formulas it is true just approximately ( Brown et al|2008 ) . In this context, the accuracy of the results 
is not affected by the use of this approximation. 



has an order of magnitude in agreement with the emission measure obtained from photon 
spectral fitting. On the contrary, the emission measure predicted by Model 3 is two orders 
of magnitude lower. 



15-Apr-2002 


Time [UT] 


EM[10 49 cm- 3 ] 


kT[keV] 


7 


A[arcsec 2 ] 


As[arcsec] 


00:03:00-00:06:00 


0.224 ±0.013 


2.03 ±0.03 


8.2 ±0.2 


388 


[-7.2,-2.0] 


00:06:00-00:09:00 


0.412 ±0.025 


1.88 ± 0.03 


8.1 ±0.2 


460 


[-7.8,-2.0] 


00:09:00-00:12:00 


0.51 ±0.04 


1.84 ±0.03 


8.3 ±0.2 


522 


[-9.2,-3.4] 


16-Apr-2002 


Time [UT] 


£M[10 49 cm- 3 ] 


kT[keV] 


7 


A[arcsec 2 ] 


As[arcsec] 


13:10:00-13:15:00 


0.331 ±0.022 


1.83 ±0.03 


8.3 ±0.1 


480 


[ 1.0, 6.2] 


13:15:00-13:20:00 


0.55 ±0.04 


1.613 ±0.021 


9.3 ±0.1 


483 


[ 1.0, 9.1] 


21-May-2004 


Time [UT] 


EM[l0 49 cm- 3 ] 


kT[keV] 


7 


A[arcsec 2 ] 


As[arcsec] 


23:47:00-23:50:00 


0.354 ±0.027 


1.85 ±0.03 


8.1 ±0.1 


215 


[ 1.0, 3.4] 


23:50:00-23:53:00 


0.62 ±0.04 


1.75 ±0.03 


8.5 ±0.1 


216 


[ 1.0, 4.4] 



Table 1: Main characteristics of the events under analysis. First column: time range con- 
sidered; second column: emission measure determined from spectroscopy; third column: 
temperature determined from spectroscopy; fourth column: photon spectral index deter- 
mined from spectroscopy; fifth column: flare area computed from the reference image at 14 



keV; sixth column: averaging interval determined as in Torre et al. (2012). 



5. Conclusion 



We have considered the influence of the background electrons drift velocity in the ac- 
celerated electrons energy loss rate due to Coulomb collisions. The background motion has 
been explained in terms of return currents that instantaneously occur in solar flares to re- 
store charge equilibrium. The obtained model for the energy loss rate has been compared 
with previous models by means of a validation process based on the use of averaged electron 
flux maps obtained from RHESSI hard X-ray imaging spectroscopy data for three different 
sources (April 15 2002, April 16 2002 and May 21 2004 ). 
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15-Apr-2002 


Time[UT] 


Mod 


E [keV] 


n[cm 3 ] 


h s [cm 5 keV 1 s x ] 


x 2 


00:03:00-00:06:00 


Mod 1 


7.7 ±0.2 


(8.0 ±0.8) x 10 10 


(1.6 ±0.3) x 10 30 


0.80 




Mod 2 


2.5 ±0.1 


(1.81 ±0.16) x 10 11 


(4.7 ±0.8) x 10 30 


1.04 




Mod 3 


7.9 ±1.9 


(7.86 ±0.12) x 10 9 


(4.2 ±0.5) x 10 30 


0.82 




Mod 4 





(3.4 ±0.3) x 10 10 


(9.52 ±1.19) x 10 30 


3.58 


00:06:00-00:09:00 


Mod 1 


8.2 ±0.3 


(6.3 ±0.6) x 10 10 


(1.04 ±0.21) x 10 30 


0.33 




Mod 2 


3.2 ±0.1 


(1.42 ±0.11) x 10 11 


(3.3 ± 0.5) x 10 30 


0.55 




Mod 3 


3.5 ±0.8 


(6.99 ±0.12) x 10 9 


(3.5 ±0.4) x 10 30 


0.43 




Mod 4 





(2.30 ±0.19) x 10 10 


(7.2 ± 0.8) x 10 30 


2.45 


00:09:00-00:12:00 


Mod 1 


4.9 ±1.2 


(1.6 ±0.3) x 10 10 


(5.1 ±2.3) x 10 28 


0.65 




Mod 2 


3.2 ±0.6 


(4.5 ±0.7) x 10 10 


(3.45 ±1.11) x 10 29 


0.70 




Mod 3 


0.2 ±0.1 


(4.47 ±0.11) x 10 9 


(1.00 ±0.24) x 10 30 


0.79 




Mod 4 





(6.2 ±0.5) x 10 9 


(1.4 ±0.3) x 10 30 


1.21 


16-Apr-2002 


Time [UT] 


Mod 


E [keV] 


n[cm 3 ] 


h s [cm 5 keV 1 s x ] 


x 2 


13:10:00-13:15:00 


Mod 1 


1.8 ±1.7 


(1.5 ±0.6) x 10 10 


(2.8 ±2.5) x 10 28 


1.22 




Mod 2 


1.2 ±0.8 


(4.9 ±1.8) x 10 10 


(2.6 ±1.8) x 10 29 


1.18 




Mod 3 


0.1 ±0.1 


(7.7 ±0.3) x 10 9 


(5.7 ±3.9) x 10 29 


1.28 




Mod 4 





(1.2 ±0.3) x 10 10 


(1.3 ±0.8) x 10 30 


1.12 


13:15:00-13:20:00 


Mod 1 


2.1 ±2.0 


(1.4 ±0.7) x 10 10 


(3.23 ± 3.24) x 10 28 


0.20 




Mod 2 


1.4 ±0.9 


(4.6 ±1.8) x 10 10 


(2.5 ±2.0) x 10 29 


0.20 




Mod 3 


0.0 ±0.1 


(7.3 ±0.6) x 10 9 


(6.8 ±4.8) x 10 30 


0.23 




Mod 4 





(9.5 ±2.3) x 10 9 


(1.0 ±0.8) x 10 30 


0.30 


21-May-2004 


Time[UT] 


Mod 


Eo[keV] 


n[cm 3 ] 


h s [cm 5 keV 1 s x ] 


x 2 


23:47:00-23:50:00 


Mod 1 


8.7 ±0.2 


(7.6 ± 0.6) x 10 10 


(1.55 ±0.23) x 10 30 


3.85 




Mod 2 


4.0 ±0.1 


(1.69 ±0.10) x 10 11 


(5.2 ± 0.6) x 10 30 


12.11 




Mod 3 


5.5 ±0.9 


(6.65 ±0.08) x 10 9 


(6.3 ±0.5) x 10 30 


10.43 




Mod 4 





(2.23 ±0.14) x 10 10 


(1.16 ±0.11) x 10 31 


40.91 


23:50:00-23:53:00 


Mod 1 


8.3 ±0.2 


(8.3 ±0.6) x 10 10 


(1.78 ± 0.24) x 10 30 


5.42 




Mod 2 


3.7 ±0.1 


(1.81 ±0.10) x 10 11 


(5.7 ±0.6) x 10 30 


17.58 




Mod 3 


3.5 ±0.5 


(7.79 ± 0.08) x 10 9 


(1.05 ±0.06) x 10 31 


21.35 




Mod 4 





(2.42 ±0.13) x 10 10 


(1.21 ±0.11) x 10 31 


54.12 



Table 2: Plasma parameters obtained with the four models described above. First column: 
data time intervals; second column: energy loss rate models used for the fit; third column: 
return currents' energies obtained from the fit; fourth column: target densities obtained from 
the fit; fifth column: averaged source amplitudes obtained from the fit; sixth column: \ 2 
values of the fit. 



11 



In all the analyzed cases, when return currents are significant to describe the observations 
(i.e., except in the case of the April 16 2002 data sets, where all models are comparable), the 
new model better fits the observed data. The main impact of this result on the theoretical 
picture of solar flares is concerned with the effectiveness of the emission process. Indeed the 
background coherent motion due to return currents tends to lower the energy loss rate of 
Coulomb collisions, and therefore the path of the electrons injected in the flare tends to be 
longer. Accelerated electrons have thus more time to emit hard X-rays by bremsstrahlung, 
which increases the efficiency of the emission process. 
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Fig. 2.— 15 April 2002, 16 April 2002 and 21 May 2004 events: R(E) as a function of 
electron energy, the blue line is the best fit hot model with return currents, the purple line 
is the cold model with return currents, the red one is the model with ohmic losses, while the 
orange line is for the hot model without return currents . 
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